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ABSTRACT 

Models for the spectra and the light curve, in the photospheric as well as 
in the late nebular phase, are used to infer the properties of the very radio- 
bright, broad-lined type lib Supernova 2003bg. Consistent fits to the light curve 
and the spectral evolution are obtained with an explosion that ejected ~ 4M 
of material with a kinetic energy of ~ 5 10 51 erg. A thin layer of hydrogen, 
comprising ~ O.O5M , is inferred to be present in the ejecta at the highest 
velocities (v > 9000km s _1 ), while a thicker helium layer, comprising w 1.25M , 
was ejected at velocities between 6500 and 9000 km s . At lower velocities, 
heavier elements are present, including ~ O.2M of 56 Ni that shape the light 
curve and the late-time nebular spectra. These values suggest that the progenitor 
star had a mass of ~ 20 — 25M (comparable to, but maybe somewhat smaller 
than that of the progenitor of the XRF/SN 2008D). The rather broad-lined early 
spectra are the result of the presence of a small amount of material (~ O.O3M ) 
at velocities > 0.1c, which carries ~ 10% of the explosion kinetic energy. No 
clear signatures of a highly aspherical explosion are detected. 

Subject headings: supernovae: general — supernovae: individual (SN2003bg) - 
stars: evolution — nucleosynthesis 
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Introduction 



Stars with an initial mass exceeding 8 — 1OM end their life when their core collapses 
to a compact remnant and the outer layers are ejected in a powerful explosion, a Supernova 
(SN). If the outer hydrogen envelope of the star is still in place at the time of the explosion, 
the SN is characterized by strong H lines and is classified as type II (see iFilippenkol 119971 . 
for a primer on SN classification). There are however cases when the outer layers of the star 
have been lost before the explosion. In particular, if the H envelope was lost but the He 
envelope was retained, the SN shows He lines and is classified as a type lb, while if both the 
H and He shells were lost neither H nor He lines are seen and the SN is classified as type 
Ic. Collectively, the class of stripped-envelope SNe is called SNelb/c. Type lb SNe are rarer 
than type Ic SNe. These subclasses should not be thought of as rigid separations, but rather 
as an indication of what is probably a continuum of pro perties. For example, evidence of 



resid ual H lines has been claimed in some SNe lb (e.g. iDeng et al.l |2000| ; lAnupama et al. 



20051 ). 



The study of stripped-envelope SNe gained momentum with the discovery t hat some ex- 
trem e ly energetic SNe Ic ( "h ypernovae" : HNe) are linked to Gamma- Ray Bursts (jGalama et al. 



19981 ; llwamoto et al.lll998l ). Thereafte r, various studies o f SNelb/c h ave been performed, 



both of individual SNe [e.g. SN 1999ex (IHamuv et a, 



20061 )] . and systematic (e.g. iMatheson et al. 



2003 



..20021); SN2004aw (ITaubenberger et al. 
Soderberg et alJbood ). " 



Some SNe Ib/c have also b een found in assoc iation with X-ray Fla shes (XRF), the softer 



and weaker vers ions of GRBs (IHeise et al.l l200lh : the SNIc 2006aj (jCampana et al.l 12006 



Pian et alJliooS ) and th e SN lb 2008D dSoderberg et aliboosh Wh i le the nature of the XRF 
is debated in both case s ( Campana et alJliood; Ghisellini et al. 2007 ; Waxman. Meszaros. fc Campana 



2007 ; Soderberg et al. 20081 ; Mazzali et al. 20081 ). both SNe turned out to be less energetic 

c than typical SNe Ib/c 
10 51 erg of energy (e.g. 



(Mazzali et al. 


2006, 


2008 




Tanaka et al. 


Iwamoto et al. 


1994 




Sauer et al. 


2006). 



It is no w clear that SNe Ic, and to so me extent also SNe lb, show a large spread of prop- 
erties (e.g. IClocchiatti &; Wheeler! 1 19971 ). Studies of various individual objec ts have begun 



to paint a p ic ture where the so-cal l ed 'typical' SNe Ib/c, su c h as SN 19941 [e.g . lNomoto et al. 



( 119931 . Il994h : IWheeler et all (ll994h ; lFilippenko et all (ll995[ ): lRichmond et al.l fll996l )] are just 
a subset of this extended family. Thus, the definition of typical may need revising. The mass 
of the progenitor stars of SNe Ib/c is not constrained, except that it must be larger than the 
minimum mass to reach core collapse (8 — 10M Q ), if the progenitor was a single star. 



The final fate of massive stars depends mostly on mass, but also on other proper- 
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ties, such as metallicity , rotation, and magnetic fields (e.g. 



Frver 


1999: 


Heeer et al. 


2003 



Woosley &: Jankall2005l ). These dependencies on often poorly known quantities add some un- 
certainty to the determination of the properties of the progenitor from those of the SN ejecta. 
Presently, SNIb/c observations and models indicate that the basic SN properties (e.g. the ki- 
netic energy liberated i n the explosion, Ek , the mass ejected, M ej -) are apparently related to 
the progenitor's mass (INomoto et al.ll2005l ). It takes very massive stars (Mzams ~ 4OM ), 



possibly of low me tallicity, to make a GRB in a highly aspherical explosion that leaves be 



hind a black hole (IFryerl Il999l ; iMacFadven fc Wooslevlll9 99). Less massive stars may still 



collapse to a black hole, po s sibly via fallback jFrverlll999h. but not give rise to a GRB [e.g. 



SNe 2002ap, Mazzali et all (120021 ) or 19 97ef. Ilwamoto et al.l fl2000h : Mazzali et al.l (1200 If )]. 



At lower masses still \M-zams ^ 22M , (lFryerlll999l ). at least if the effect of rotation is not 



considered], the star collapses to a neutron star and may eject some material at relativistic 
velocities, but in this case t he relativistic outflo w is less energetic and is observed as an 

The se less energetic explosions may be 



X-Ray Flash [e.g SN2006ai (IMazzali et al.l 120 06) 
less asymmetric ( Mazzali et alj|2007l ; iMaeda et al. 



20081 ). 



One general feature of SNe Ib/c is that their progenitors must have lost at least the 
outer hydrogen envelope. If a hydrogen envelope more massive than ~ 0.1M Q is left, the 
star would probably have a RSG structure at the time of core collapse, and the SN light 
curve would be influenced initially by hydrogen recombination and appear as a type II light 
curve. Some transitional objects like SN 1993J only had a very small hydrogen envelope, 
whose presence did not affect the light cu rve significantly but was responsible for the Type 
II spectrum this SN displayed early on (INomoto et al.lll993l ). Spectroscopically, after an 
early Type II phase SN 1993 J turned to a Type lb, and is therefore known as an SNIIb. 
Only few SNe lib have been observed, and little is known about their spread in mass and 
energy. 

Exactly how a star manages to preserve just a thin shell of H is debated. The scenarios 
are different depending on whether the star was single or a member of an intera cting binary. 



In th e former case, a mass lo ss occurs via a radiatively driven stellar wind (e.g. IWeiler et al. 



20071 ; IStockdale et al.l 120071 ) . This can be effective at relatively high masses (> 30M Q ). 
Stars with masses of ~ 27 — 28M Q may end their lives with only a thin H s hell left in place 



dCrockett et al.l 120081 ) . Such stars would be classified as WRL or WRH (ISmith &: Conti 
20081 ). However, this mechanism is much less efficient for less luminous stars. Yet, the infer- 
ence for events such as SN 1993J is that the progenitor had a mass of ~ 15M . In this cas e 



the la tter scenario, binary interaction, seems more appealing. For example, INomoto et al. 



(119951 ) suggest that an SNIIb is a result of the spiral-in of a low-mass companion into a 
massive star, which causes the ejection of most of the massive star's H-rich envelope. This is 
probably a good scenario to produce an energetic explosion, because it can explain both how 
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the core gains angular momentum from spiral-in, and how most of the envelope is ejected 
by spiral-in heating. There is some observational evidence in support of the bina ry scenario: 



for SN 1993J a spectroscopic companion has been observed (IMaund et al.ll2004f). Ev idence 



of a binary companio n has also been c laime d for the SNIIb 2001ig (IRyder et al 



20061 ). while 



for the SN lb 2008ax ICrockett et al.l (120081 ) suggest that binary interaction may be a more 
viable scenario th an that of a single star. A similar debate exists also regarding the origin 
of SNe Ib/c (e.g. iNomoto et al1ll993l . boOlh . 



SN2003bg was initially classified as a broad-lined SNIc (jFilippenko Chornockll2003l ). 
but it soon developed a broad Ha P-Cygni feature as well as He I lines, making it a SNIIb. 
This is the first SNIIb for which broad lines have been observed, suggesting that it is 
a SNIIb hypernova . The optical data of SN2003bg are presented in a companion paper 
( iHamuy et al.ll2009l ). In the nebular phase the spectra were similar to those of SNe Ib/c, but 
evidence of H may be present in the form of a weak, flat-topped Ha emissio n, reminiscent of 



SN 1 993J. Additionally, SN2003bg was exceptional for its radio luminosity. ISoderberg et al. 



( 120061 ) suggest that the observed radio light curve can be interpreted as the signature of 
sub-relativistic ejecta (with velocity v ~ 0.2c) producing synchrotron radiation as they im- 
pacted into a medium that may have had the properties of a WR-star wind (mass loss rate 
M ~ 4 lO _5 M yr~ 1 , wind velocity v ~ 1000km s -1 ). They also suggest that the wind 
may have been highly variable, especially just before the explosion. This is reminiscent of 
the situation of t he XRF-SN2006ai, which was probably s urroun ded by a thick wind shell 



Campana et al.l (120061 ); IWaxman. Meszaros. fe Campanal (120071 ). but see iGhisellini et al. 



( 120071 )] . Interestingly, the radio light curve of SN2003bg and that of SN2001ig were very 
similar. 

The available evidence therefore suggests some similarity among all SNe lib, despite 
clear differenc e s in t he ejected mass and the kinetic energy of the explosion: for SN 1993 J 
Nomoto et al.l (119931) estimated an e j ected mass of 4M and an explosion E K of 10 51 erg, 
while for SN 2001ig ISilverman et al.l ( 20091 ) estimate an ejected mass of ~ 1M Q below a 
velocity of 4300 km s _1 . 

Here we concentrate on extracting the basic properties of SN 2003bg, and hence of the 
progenitor star and of the explosion. In order to do this we model the light curve and the 
spectra in both the early photospheric phase and the late nebular phase. 
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2. Data 

Observations of SN 2003bg started very soon after the putative explosion date (2003 Feb 
25 UT), an d cover the evolutio n of the SN for almost one year, extending into the nebular 



phase (see lHamuy et al.l 120091 . for more details). The peak phase was unfortunately not 
covered, and there is a gap in the data at the onset of the linear decline phase. Still, the 
available data make it possible to constrain the properties of the SN reasonably accurately, 
especially in view of the availability of both early and late spectra. 

We have modelled 10 spectra, of which 8 were obtained around the peak and 2 in 
the nebular phase. Combining the analysis of the outer part of the ejecta conducted by 
means of the early-time spectra and the determination of the properties of the inner, slow- 
moving ejecta made possible by the nebular spectra, both the density and the abundance 
stratification in the SN ejecta can be determined. 



3. Models for the early phase 

Unlike Type la SNe, the modelling of SNe Ib/c must establish both Ek and M e j simul- 
taneously. As these are the two parameters that play the major role in determining both 
the light curve and the spectra, clearly a reliable solution must be able to reproduce simul- 
taneously both the spectral evolution and the light curve. As is well known, the solution o f 
the light curve alone is degenerate to Eg and M e j if both are allowed to vary (jArnettlll982l ). 
as is the case for SNe Ib/c. 

The first steps to obtaining a solution are therefore to consider that: i) the light curve 
has a rather narrow peak: the presence of H and He does not influence the light curve much, 
because the H shell must be thin and He does not contribute significantly to the opacity. 



(Mazzali et al. 


2002) 


~ O.2M of 56 Ni. 1 
trypernovae such as 


(Mazzali et al. 


2006) 



2008 



Tanaka et all 120091 ) . 



3.1. Method 



The early-time spectra were modelled usi ng a Monte Carlo radiative transfer code 
f lMazzali fc Lucvlll993l : lLucvlll999l : lMazzah1l2000f ). The code assumes that the SN luminosity 
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emerges from a sharp, grey photosphere (Schuster-Schwarzschild approximation). This ap- 
proximation is valid if most 7-rays and positrons produced in the radiactive decay of 56 Ni to 
56 Co and hence to 56 Fe deposit their energy, which is eventually thermalised and converted 
into optical photons, below the photosphere. This is the case in the early phases of a Type 
I SN. Using as input the SN luminosity and the velocity of the photosphere at the desired 
epoch of the model (velocity and time give the photospheric radius since the expansion of SN 
ejecta is homologous: v = rt), an approximate effective temperature at the photosphere can 
be obtained assuming black body emission (L = 4tiR 2 (JbT 4 ', where Ob is Stefan-Boltzmann's 
constant). 

Photons emitted at the photosphere propagate in the expanding ejecta, where radiative 
equilibrium is assumed. Photons can undergo electron scattering and they can in teract 
with spectral lines. In the latter case, photon branching is explicitly considered (ILucy 
19991 ; iMazzalil |2000| ). This process allows photons to shift from wavelengths where line 
opacity is very high to spectral regions where line opacity is low, and thus to escape the 
SN envelope, and it is essential in shaping the spectra of Type I SNe. Level populations 
and the ionization structure in the SN envelope are computed taking into account deviations 
fr om Local Therm o dynam ic E quilibrium (LTE ) using a nebular approximation as discussed 



111 



Abbott fc Lucyl (119851 ) and iMazzali fc Luc y (1993 ). This appr oximation has been shown 
to reproduce well the effects of non-LTE (IPauldrach et al.lll996l ). Iterating a series of MC 



experiments, a temperature structure in the ejecta is then c omputed. F inally, the emerging 
flux is recorded using an efficient formal integral approach (ILucyi Il999h . The calculation of 
the observed flux requires the values of the distanc e and the reddening to the SN. The code 
has been applied to a number of Type I SNe (e.g. IMazzali et al.ll2008l ). 



Apart from the parameters described above, and for atomic models, the only other 
input required by the code is the run of density in the ejecta with expansion velocity. This is 
known as an explosion model. The properties of the model determine the appearance of the 
spectra. The basic quantities that define an explosion model are the ejected mass and the 
kinetic energy of the explosion. In particular, the presence o f large amounts of material at 
high velocities leads to spectra characterized by broad lines (IMazzali. Iwamoto. &: Nomoto 
2000 ]), while, given the same kinetic energy, more massive models have broader light curves 
The simultaneous fit of the light curve and the spectra of an SN using one 



explosion model therefore constitutes a validation of the model itself. 

Because of the spectral similarity with SNe 1997ef and 2002ap we started with an 
explosion model based on CO100, a model characterized by M e j~ 8M and E K ~ 10 52 erg. 
Model CO 100 was used to reproduce the light curve and the spectral evolution of SN 1997ef, 
but the very broad lines observed in this SN required the addition of a high-velocity, low- 
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density tail which increased the Ek to ~ 2 10 52 erg. The model was also used in rescaled 
form (1/4 of the mass) to fit the SNIc 200 2ap, which was char acterized by broad spectral 
features but a relatively narrow light curve (IMazzali et al.ll2002l ). SN2003bg poses a similar 
problem. In order to reproduce its narrow light curve we needed to rescale model CO 100 
down in mass. Furthermore, line blending in SN2003bg is less extreme than in SNe 1997ef or 
2002ap - it is more similar to SN2006aj: e.g. the Ca II IR triplet and the O I 7774 A line do 
not blend - and so the density gradient in the outer part must have been steeper. Therefore 
we did not need the high- velocity tail, and we finally used a model with half the mass of the 
original CO100. The model we selected has M e j ~ 4M Q and Ek ~ 5 10 51 erg. 

In the spectral models the luminosity L, which is a required input, was determined such 
as to match the flux in the spectra. Another required input is the photospheric velocity 
v p h- This was chosen so that the position of the absorption lines and the overall spectral 
distribution in the synthetic spectra reproduced the observations as closely as possible. The 
abundances in the ejecta, also a required input, were chosen to optimize the fit, but our aim 
was not so much to reproduce the details of the spectra but rather to explain the general 
evolution of the SN. 

Because H and He lines appear in the spectra we assumed that these elements are 
present in the outer part of the ej ecta. He lin es are notoriously not visible unless the levels 
of He I are non-thermally excited (jLucyi 1 1 99 lh . Non-thermal excitation is caused by the fast 
electrons which are produced in the thermalization process of the 7-rays that are emitted 
by the radioactive decay of 56 Ni into 56 Co and hence into 56 Fe, which is the source of light 
for SNe. At early times, the density in the ejecta is high and thermalization occurs near the 
place where the 7-rays are emitted. Since the helium layer is located at a larger radius than 
the 56 Ni non-thermal processes do not occur. Only at later times do fast electron begin to 
travel considerable distan ces before being ther malized. In fact, He I lines develop rather late 
in SNe lb (e.g. SN 2008D. IMazzali et alJl2008h . 



We do not explicitly treat the propagation of the fast electrons, which also requires a 
detailed mapping of the distribution of 56 Ni in what may have been an aspherical explosion. 
Therefore, the departure coefficients caused by non-thermal processes with respect to the 
level population computed with our pseudo-NLTE treatment are not explicitly computed. 
Instead, we mimic the effect of overexcitation and ionization caused by non-thermal processes 
using an arbitrary 'non-thermal factor' /hpt to multip ly all He I excited level populations. 
These factors increase with time (IMazzali fc Lucyill998l ) . reflecting the increasing penetration 
of fas t elect rons in the SN envelope, and they were taken to be similar to those derived by 
Lucy! (Il99ll ). While this is a rough treatment, it g ives one at least a feeling for the presence 
or lack of He in the ejecta (jTominaga et al.l 120051 ). A non-thermal factor was also used for 
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hydrogen (/h): since the H mass must be small, this was required to match the observed 
line strengths. 

The explosion date was chosen to be 2003 February 25. This gives a rise time to 
maximum of ~ 15 — 20 days, similar to SN 1998bw. We used a distance modulus of 31.68 
(d = 24Mpc), and a very small reddening [E(B — V) = 0.02]. The early-time spectral fits 
are presented in Figures 1 and 2 and are briefly discussed in the following. The input values 
for the spectral models are summarized in Table 1. 

3.2. Results 

3.2.1. 2003 Feb. 28 

The model is shown in Figure la. It was computed using logL = 41.47 [erg s -1 ], 
v p h— 22000 km s™ 1 , and t = 3 days. The mass above the photosphere is ~ 0.2M Q , and 
the radiation temperature at the photosphere is ~ 7000 K. Although He dominates the 
composition (90% by mass above the photospheric velocity), and a non-thermal factor / He i = 
10 3 was introduced, no strong He lines are visible, and so the He distribution cannot be 
constrained. Hydrogen is present only in a very small amount (4% by mass). If more H is 
included, the increased electron density changes the ionization balance of all other species. 
Therefore, the same non-thermal excitation factor as for He I was used also for H I (/h = 10 3 ). 
Both H and He lines make only a minor contribution to the troughs near 6000 and 4300A, 
which are dominated by Si II and Co II lines, respectively. Other features are due to Fe II, 
Si II, O I, Ca II, and Mg II and are marked in Figure la. The presence of a small amount 
of material at high velocity (~ 0.03M Q ) are located at v > 0.1c is sufficient to give rise to 
broad lines. 



3.2.2. 2003 March 2 

The model shown in Figure lb was computed with logL = 41.82 [erg s _1 ], v p h= 
18000 km s _1 , and t = 5 days. It is similar to the previous one, but the lines are deeper as 
the photosphere moved rapidly inwards. The mass above the photosphere is now « O.4M . 
The abundances are also similar, with He the dominant species. The higher L leads to a 
higher radiation temperature near the photosphere, T ~ 7500 K. Hydrogen lines now make 
an important contribution: in the feature near 6200 A, Ha is as strong as Si II 6355A, 
which contributes to making that feature broad. The same non-thermal factors as for Feb 
28 (fuel = 10 3 , /h = 10 3 ) were used for both H I and He I. He I lines are present in the 
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spectrum. He I 5686 A is blended with Na I D, which is however weaker. The inconsistent 
strength of Ha, which is too strong, and H/3, which is comparable to the observation, may 
suggest that a full, rather than a parametrized treatment of NLTE and non-thermal effects 
for at least H may be required. If the contribution of Si II 6355A to the Ha-Si II feature is 
reduced, the synthetic line does not match the velocity of the observed absorption. 

3.2.3. 2003 March 4 

This model is shown in Figure lc. It was computed with logL = 41.95 [erg s _1 ], 
v p h= 14100 km s _1 , and t = 7 days. It has a similar composition as the previous model, and 
the same non-thermal factors for H I and He I (/hcI = 10 3 , /h = 10 3 ). As the temperature 
continues to increase (now T ~ 8000 K) with the luminosity, and more mass (now ph O.7M ) 
is located above the photosphere, the H and He lines make an increasingly important con- 
tribution. Ha is now stronger than Si II 6355 A, H/9 is visible, and so are some of the He I 
lines (especially 4471 A and 5876 A). Apart from this, there are no major differences with 
respect to the previous spectrum. 

3.2.4. 2003 March 12 

This is the spectrum nearest the time of maximum, the epoch of which was not covered. 
It was computed with logL = 42.31 [erg s -1 ], v ph — 8800km s _1 , and t = 15 days, and it is 
shown in Figure Id. The spectrum is now bluer than before, but the lines are the same. The 
temperature has further increased (T ~ 8600 K), and now H I and He I lines dominate (the 
composition is similar to the previous epochs, despite the much larger mass located above 
the photosphere, ~ 1.5M , and the same non-thermal factors were used for both H I and 
He I: fuel = 10 3 , fn = 10 3 ), while Fe makes a smaller contribution by comparison. 

3.2.5. 2003 March 31 

This spectrum was obtained after an observational gap near the time of maximum. 
It is much redder than all pre-maximum spectra, which is the result of the decrease in 
temperature (T ~ 6000 K) and luminosity. The model was computed with logL = 42.29 
[erg s _1 ], Vph— 6000 km s _1 , and t = 34 days, and it is shown in Figure 2a. The mass above 
the photosphere is now 2.4M Q . The abundances are not very different, with a slight 
decrease of He, which is reduced to ~ 70% and replaced mostly by Ne near the photospheric 
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layers. The abundance of Fe-group elements has now increased to ~ 4%, which is required 
to block the near-UV flux. 

Several He I lines are now clearly visible, emphasizing the lib nature of the SN. They can 
be reproduced using / He i = 10 6 . An increase with time of the departure c oefficients can be 



expec ted because of the increased penetration of the 7-rays, as discussed by lMazzali fc Lucy 



(1l998l ). Although any statement about the He mass is made highly uncertain by the presence 
of the huge non-thermal correction, the strength of the He lines suggests that the He mass 
should be large. In any case, the He- shell is limi ted to ~ 2M at most stellar masses in 



evolutionary models (INomoto fe Hashimoto Ill988l ) 



The H I lines are now narrower, and have a smaller blueshift than lines of other elements. 
In order to get a reasonable match to these lines (only Ha and H/5 are seen) it must be 
assumed that H is present only at v > 9000 km s _1 . The non-thermal coefficient is only 
marginally larger than at previous epochs (/h = 2 10 3 ). Given the low abundance of H (less 
than 4%), we derive a total H mass of ~ O.O5M . The mass of H cannot be much larger for 
two reasons. Firstly, if the mass was much larger than 0.1M Q the star would not be compact 
and one would expect a type II light curve, with a phase dominated by H recombination. 
Secondly, a large H mass would imply a higher electron density, which would then result in a 
different overall ionization regime and in a worse fit to the spectra. The oxygen abundance 
near the photosphere cannot be determined because the spectra do not cover the O I 7774 A 
line, so the same value as in the previous epoch was used (3%). 



3.2.6. 2003 April 4 

This spectrum was obtained just a few days later and is similar to the previous one. 
The model (Figure 2b) was computed with logL = 42.24 [erg s _1 ], v p h= 5600 km s _1 , and 
t = 38 days. The mass above the photosphere is now ss 2.5M Q , and the temperature has 
decreased only slightly, to 5600 K. The lines are now more clearly separated, with metal 
lines making a larger contribution to the spectrum, although Ha, H/3, and some He lines 
are still very strong. The abundances have not changed much. The non-thermal co e fficien t 



for He is larger (fuel = 10 7 ), but comparable to what was used in iTominaga et al.l (120051 ). 
The non-thermal coefficient for H is the same as in the previous epoch (/h = 2 10 3 ). H is 
present only above 9000 km s _1 , but now a confinement for He has also been introduced at 
6500 km s _1 . This is required to fit the shape and position of the He I lines. 
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3.2.7. 2003 April 9 

The next spectrum is similar to the previous one. The model (Figure 2c) was computed 
with logL = 42.21 [erg s^ 1 ], v p h= 5100km s _1 , and t = 43 days, and has a mass above 
the photosphere of ~ 2.65M . The spectrum is marginally cooler than the previous one 
(T « 5300 K, but the lines that are present are the same. H I and He I lines are sharp. 
The abundances have changed consibrably with respect to the previous spectrum. The 
abundance of He is down to ~ 50% by mass, and O, Ne, and Si are highly abundant near 
the photosphere. The confinements of H (above 9000 km s _1 ) and He (above 6500 km s _1 ) 
are as in the previous spectrum, but the non-thermal factors are higher: /fiei — 10 9 and 
/h = 3 10 3 . Dep arture coefficients that increase with time are expected as 7-rays propagate 



more efficiently (IMazzali fc Lucylll998l ). but the precise values should be tested with realistic 
calculations. These require a detailed knowledge of the distribution of 56 Ni as well as of its 
mass, and are postponed to future work. Still, a slower rise of the non-thermal factor for H I 
than for He I may be justified by the fact that hydrogen is located further out than helium 
in the ejecta, and hence it is more separated from the source of fast electrons. Also, H I 
levels are separated by smaller energies than those of He I, and thus are more coupled to the 
thermal pool. This also applies to the ionization potentials of the two ions. 



3.2.8. 2003 April 10 

This is the last of the photospheric epoch spectra. The model (Figure 2d) was computed 
with logL = 42.20 [erg s" 1 ], v p h= 5000 km s _1 , and t = 44 days. The composition, non- 
thermal factor and confinements are the same as in the previous spectrum, which was only 
one day earlier. The mass above the photosphere is now rj 2.7 M g , and the temperature is 
~ 5300 K. H land He I lines are very sharp. The volume sampled by last few spectra in the 
photospheric epoch is just outside the volume sampled by the nebular spectra. 



3.3. Results of early-phase modelling 

The main result of the photospheric-epoch models are that SN 2003bg was a very ener- 
getic explosion (a hypernova) and that H and He are both present, but are clearly confined 
in velocity: H above 9000 km s _1 and He above 6500 km s _1 . This explains the weakness of 
the H and He lines in the early spectra and the almost complete absence of H in the nebular 
spectra: these elements are removed from 56 Ni, the source of the fast particles responsible 
for non-thermal processes. Thus the SN looked like an SNIc early on, although H and He 
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dominate the composition of the outer ejecta. This is because the H content should be small 
(only ~ 0.05M Q ) and the He I levels giving rise to the strongest lines are not excited early 
on. These levels are excited mostly via non-thermal processes, but the fast particles that 
are responsible for the excitation are produced deeper in the ejecta, and cannot penetrate 
from the deeper regions when the densities are still too hig h. The delayed develop ment of He 
l ines has been observed before in SNelb [e g. SN2005bf JAnupama et aD l2005h: SN2008D 
rtSoderberg et alihoosl : llMazzali et aD I2OO8I : kanaka et al.ll2009h . Below 6500 km s _1 , silicon 
and oxygen dominate the composition, and Fe-group elements are also significant. 

This behaviour of the H I and He I lines shows that if H and He are present in the ejecta 
they will at som e point manifest themselv es in the spectra. Since in SNelc He lines are 
never seen (e.g. iTaubenberger et al.ll2006l ). and H lines are not seen in most SNelb, these 



two subtipes must contain very little He and H, respectively. 



3.4. Models for the nebular phase 



In order to assess the properties of the inner parts of the ejecta, and thus to determine 
both the 56 Ni mass (in addition to the light curve results) and the ejected mass, it is essential 
that late-phase spectra are modelled as well as early-time ones. At advanced epochs, the SN 
ejecta become optically thin and behave like a nebula. The gas is heated by the deposition 
of the 7-rays and positrons that are emitted in the decay of 56 Co to 56 Fe. Collisions of 
the fast particles that are produced in the thermalization process excite the gas, which is 
cooled via mostly forbidden-line emission. Th e gas conditions are modelled with a non 



local thermodynam ic equilibrium (NLTE) code (lAxelrodlll980l ; iRuiz-Lapuente fc Lucylll992 
Mazzali et aPboOlh . 



The late-time spectra of SNelb/c are dominated by a strong O I 6300, 6363 A emission 
line, the shape of which can be used to inf er the geometry of t he ejecta and our orientation, 
thanks to the transparency of the nebula (IMazzali et al.ll2005l ). 



We modelled the two latest spectra of SN2003bg, obtained on 2003 Dec 16 and 23, 
corresponding to fiducial epochs of 293 and 300 days, respectively. The two spectra are 
rather similar. Besides O I 6300, 6363 A, other strong lines are Mg i] 4570A, several [Fe n] 
lines which testify to the brightness of the SN, and Ca II lines. 

The shape of the emission lines does not suggest the presence of strong asymmetries, 
and the spectra can be modelled assuming a single density and homogeneous abundances 
below a velocity of 5000 km s _1 , as indicated by the width of the strongest emission lines 
([O i], Mg i]). Models for the two spectra yield reasonably consistent results (Fig. 3). 
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Given the somewhat uncertain flux calibration we take the mean of the results of the two 
spectral fits. The spectra can be reproduced assuming an emitting mass of ~ 1.6M located 
below a velocity of 5000 km s _1 . The main constituent of the emitting nebula is oxygen 
(~ O.9M ). Other elements are C (« 0.1M o ), Si (« 0.3M Q ), S (« O.1M ), Ca (» O.O2M ), 
Mg (ps O.OO5M ), and 56 Ni. The mass of 56 Ni in the nebula is ps O.16M , as required by the 
simultaneous fit of the [Fe n] and other lines. By the time the nebular spectra were observed 
most 56 Ni had decayed to 56 Co and 56 Fe. 

The explosion model we used for the early-time models has a mass of 1.3M below 
5000 km s -1 . The result from the nebular spectra is slightly larger, indicating the presence 
of more mass at low velocities than included in a spherically symmetric ex plosion model. 



This i s often seen in SNe lb/c [e.g. SN 19941 flSauer et al.l 120061 ): SN 1998bw flMazzali et al. 



2001 



Maeda et al. 1 IgOOgl )]. and may be the result of some asphericity in the explosion 



flMaeda fe Nomotol 12003^. which is typ ical of HNe and possibly of all core-collapse SNe 



^Leonard et alJbood l Maeda et al.l 120081 ) . Indeed, a hint for the presence of an inner core of 



material with higher density may be found in the profile of the [O i] 6300, 6363 A line, which 
may be composed of a main profile characterized by nar rowly separated dou ble peaks plus 
a central core and is reminiscent of that of SN2002ap (IMazzali et al.l 120071 ). A model in- 
cluding a density discontinuity to reproduce this profile in detail may have a slightly smaller 
total mass. Still, the fact that the deviation from parabolic profiles is very small in all lines 
suggests that any asphericity is small and probably affects only the innermost ejecta. 

The strength of the semi-forbidden Ca n] 7200 A line with respect to the allowed Ca II 
IR triplet suggests that the density was low. In our models, the electron density is n e ~ 
2 10 6 cm -3 . A weak, broad [O i] 5557 A line may be present in the data, and is only partially 
reproduced in the synthetic spectra, which also suggests a slightly higher density. A weak, 
broad emission to the red of [O i] 6300, 6363A may be identified as Ha, as in SN 1993 J 
jPatat. Chugai. fe Mazzalilll995h . 



3.5. Results of spectral modelling 

Combining the early- and late-time modelling results we find that a model of a highly 
energetic explosion which ejected ~ 4 — 5M of material with EkPz 5 10 51 erg can reproduce 
the observations. Only a small mass of H is present, ~ O.O5M , located at v > 9000 km s _1 . 
A massive He shell is located at v > 6500 km s _1 , comprising ~ 1— 2M . The ejected CO core 
mass is ~ 3M , of which ~ 1.3M is oxygen. The rest is mostly carbon (~ O.15M ), neon 
(~ O.5M ), silicon (~ O.6M ), and sulphur (~ O.2M ). The mass of 56 Ni is ~ O.17M , most 
of which is located inside of 5000 km s _1 . The confinement of 56 Ni to the lowest velocities 



-14- 



explains the late excitation of the He lines. 



Layers of H and He at the highest ejecta velocities were found by iBranch et al.l (120021 ) 
in most SNe lb. In SN 2003bg and in other SNe lb the H mass must be larger, so that it gives 
rise to strong spectral lines. Also, in the case of SN2 003bK, we find that H extends d eeper 
than in most SNe lb, similar to the type lib SN 1993 J (IPatat. Chueai. fc Mazzalilll995l) . The 
inner distribution of helium is at the lower limit of the Branch et al. ( 20021 ) sample. 



Early on, neither H no r He lines are clearl y detected. An elegant explanation of this 
behaviour of H was given by lBranch et al.l (120021 ) . Additionally, in the case of SN 2003bg, the 
early identification of H lines is rendered difficult by the breadth of the spec tral features. He 
lines only develop later, as non-thermal excitation increases with time (see iMazzali fc Lucy 
1998h . 



The velocity evolution of SN2003bg (Figure 4) is similar to that of other SNelb/c, in 
particular the HN SN2002ap, although our estimate of the explosion epoch is uncertain. 
It is also similar to SN2008D early on. The velocities we derived are based on modelling 
the spectra, and are therefore influenced by trying to reproduce both the overall spectral 
distribution and the observed line width. Therefore, they are sensitive to the assumed epoch, 
especially at the earliest phases. This is an indirect way of assessing the reliability of our 
choice of epoch. If the date of explosion was chosen to be earlier than assumed here, the 
velocity required to fit the spectra would most likely be smaller, and the synthetic lines 
would be less blended. A later epoch of explosion, on the other hand, would result in even 
larger velocities, especially at the earliest times, and the blueshift of the lines would be too 
large. 



4. Light Curve Model 

We synthesized the bolometric light curve (LC) of the ejecta models as constrained by 
spectral fitt ing, and compare t hem to the observed LC, which is described in the accompany- 



ing paper (iHamuy et al.ll2009l ). We used the 1-D SN LC code that was originally developed 
by Iwamoto et al. (2000). The code solves the energy and momentum equations of the radia- 
tion plus gas in the co-moving frame, and is accurate to first order in v/c. Electron densities 
and the electron scattering opacity are determined from the Saha-Boltzmann equation. The 
energy deposition from radioactive decays of the newly synthesized 56 Ni and its daughter 
56 Co is calculated with a gray 7-ray transfer code, assuming an absorptive opacity of 0.05Ye 
cm 2 g _1 (Swartz, Sutherland & Harkness 1995). Unlike the sophisticated but time-consuming 
recipes in the original code, we adopted the approximation proposed by Gomez-Gomar & 
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Isern (1996) for the Eddington factors, and fitted the TOPS opacities (Magee et al. 1995) 
to find an empirical relationship between the Rosseland mean and the electron scattering 
opacity. Such simplifications were also made by the authors when modeling the LCs of other 
Type Ic SNe. 

The best-fitting model LC shown in Fig. 5 (solid line) is characterized by an ejected 
mass of ~ 4.8M , with a kinetic energy of ~ 5 10 51 erg and including O.15M of 56 Ni, whose 
distribution in the ejecta was adjusted by hand. The density structure above 5000 km s _1 was 
obtained scaling model CO100 as constrained by our photospheric-epoch spectral models. 
Below 5000 km s _1 , the density was enhanced to the level required by the nebular spectrum 
modelling. In order to retain the calculation accuracy, the ejecta material with v/c > 0.1 was 
cut out. This amounts to a negligible mass because of the extremely low density at such high 
velocities. Except for 56 Ni, we adopted the abundance distribution of other elements derived 
from the spectrum synthesis. In particular, H exists only at v > 9000 km s _1 . The total 
H mass is only 0.05 M , which is too small to affect the LC shape by H recombination, 
= as expected. The model LC shown has a 56 Ni distribution of ~ 0.01M Q (or 1% in 
mass fraction) between 10,000 and 20,000km s^ 1 , « 0.04M Q (« 2.5%) between 5,000 and 
10,000 km s _1 , « O.O5M (« 4.4%) between 3,000 and 5,000 km s~\ and « 0.05M Q (« 4.4%) 
below 3,000 km s _1 . Thus the total 56 Ni mass derived from the light curve (0.15M Q is similar 
to that obtained from the spectral fits, but the distribution is different. The light curve model 
requires more 56 Ni at high velocities. 

For comparison, two other model LCs with different model parameters are also shown 
in Fig. 5. The dashed line shows a LC obtained with the same density structure and hence 
the same mass and kinetic energy as the best-fit LC, but with the 56 Ni distribution derived 
from the spectrum modelling. The total 56 Ni mass is then « O.18M , of which ps 0.16M o 
are located below 5000 km s _1 . This low-velocity 56 Ni must be responsible for the high peak 
and post-peak flux, which is above the early-time photometry by ~ 0.3 mag, although the 
LC matches the very late photometry, as it should by design, since the inner distribution of 
56 Ni was obtained from fitting the nebular spectra. Moreover, there are only O.OO2M of 
56 Ni above 10,000 km s _1 , not enough to power a very fast early LC rise. 

The dotted line is the synthesized LC of a model without the density enhancement 
below 5000 km s _1 and with the 56 Ni distribution constrained by spectral modeling. The 
total ejecta mass, kinetic energy, and 56 Ni mass are ~ 3.9M , 5 10 51 erg, and ~ O.12M , 
respectively. Although the energy available is smaller, the LC peak is actually as as bright 
as in the other models because less ejecta material is heated. However, because of the small 
ejecta mass relative to the kinetic energy, the LC peak is too narrow and the LC drops too 
rapidly when compared to the observations. As a further test of our model, we plot in the 
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inset of Fig. 5 the evolution of the photospheric velocity of the three models as estimated 
approximately by the LC code. They all reproduce well the values that were determined by 
spectral modeling (stars). 

While we cannot exclude the possibility that other 56 Ni distributions, with some tuning, 
may also give a good LC fit, it would be a futile attempt to try to pin down the 56 Ni dis- 
tribution so precisely in a one-dimensional model. Nonetheless, the distribution we adopted 
results in a rapid LC rise, and in a maximum luminosity and a broad LC peak that match 
the observations well. After about day 150 the model LC lies below the observations by 
~ 0.2 mag, but the slope is comparable. 



5. Discussion 



The properties of SN2003bg derived from our analysis suggest that the SN was the 
explosion of a massive star, clearly more massive than the Type lib SN 1993J. The large 
explosion energy (Ek ~ 5 10 51 erg) qualifies SN 2003bg as a hypernova, in line with the broad 
lines detected in the earliest spectra. If the ejected mass was ~ 4 — 5M Q , of which ~ 1 — 2M 
were He and only ~ 0.05M© was H, the star had lost only its H envelope. Assuming that 
the compact object formed at the time of core-collapse had a mass between 1.4M (if it 
was a neutron star) and 2M (if it was a black hole), we have a He star of ~ 5.5 — 7M & 
at explosion. This would imply a CO core of ~ 4 — 5M and a star of main- sequence mass 
Mzams ~ 22 — 25M , very similar to the star that exploded as SN2002ap (IMazzali et al. 



20061 ) and also similar to, or perhaps somewhat less massive than the progenitor of the 
XRF-SN 2008D (IMazzali et al.ll2008l : kanaka et al]l2009h . SN 2003bg had a similar explosion 



energy as these two other SNe, but because of the larger M e j very high velocities were not 
reached, and spectral lines were not as broad. This situation (low velocities, larger overlying 
mass) also would have made the generation of a hard X-ray or 7-ray transient very unlikely. 

The 56 Ni distribution derived from the light curve model is different from that obtained 
from the spectra. Most of the 56 Ni obtained from the spectral models is required to power 
the nebular spectra. If a higher abundance of 56 Ni was used outside of the 5000 km s _1 zone 
sampled by the nebular spectra, the synthetic [O 1] 6300A line would be too broad, unless 
56 Ni and oxygen were physically separated. This would be a sig nature of asymmetry in the 
explosion, as was seen in, e.g. , SN 1998bw (IMazzali et al.ll200ll ). 



SN 2003bg falls on the relations between stellar main se quence mass and, re spectively, 
mass of 56 Ni synthesised and Ek produced in the explosion (jNomoto et all 120051 ). as shown 
in Figures 6 and 7, suggesting that these relations hold not only for SNe Ic but also for less 
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stripped core-collapse SNe. 

Given our estimate of the progenitor mass, the compact remnant is likely to have been 
a black hole, but it also could have have been a neutron st ar. Overenergetic explos ions are 



expected also in this case if the neutron star is a magnetar (jBucciantini et all 120081 ) . 



However, the pres ence of significant CSM gave rise to strong radio emission in SN 2003bg 
JSoderberg et alibood ). What is the difference between SN 2003bg and SN 2006aj? SN 2006aj 



had lost its He envelope, and possibly part of its CO core. Thus the ci rcumstellar materia l 



inside of which SN2006aj exploded was composed mostly of oxygen (IMazzali et all 120061 ) 



SN2003bg only lost m ost of its H envelope, and H was the main constituent of the CSM 



(ISoderberg et al.ll2006l ). Quite possibly, SN2006aj was a binary star, which helped ejecting 
the envelope. SN 2003bg may have been a single star, and if it was an interacting binary the 
effect of interaction was not as dramatic as in SN2006aj. 
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Table 1. Parameters of the early-time synthetic spectra 



Date 


SN epoch 


logL 


Vph 


/(He I) 


/(H I) 




[days] a 


[erg s _1 ] 


km s _1 






28 Feb 2003 


3 


41.47 


22000 


10 3 


10 3 


2 Mar 2003 


5 


41.82 


18000 


10 3 


10 3 


4 Mar 2003 


7 


41.95 


14100 


10 3 


10 3 


12 Mar 2003 


15 


42.31 


8800 


10 3 


10 3 


31 Mar 2003 


34 


42.29 


6000 


10 6 


210 3 


4 Apr 2003 


38 


42.24 


5600 


10 7 


210 3 


9 Apr 2003 


43 


42.21 


5100 


10 9 


310 3 


10 Apr 2003 


44 


42.20 


5000 


10 9 


310 3 



a Thc epoch is given from the putative date of explosion, 25 Mar 2003. 




Fig. 1. — Observed spectra of SN 2003bg (green/black) compared to synthetic models 
(black/grey) for pre-maximum epochs. 
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Fig. 2. — Observed spectra of SN 2003bg (green/black) compared to synthetic models 
(black/grey) for post-maximum epochs. 
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Fig. 3. — Observed nebular spectra of SN 2003bg (black) compared to synthetic models 
(red/grey). 
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Fig. 4. — The evolution of the photospheric velocity obtained from spectral models in a 
number of SNe Ib/c. 
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Fig. 5. — The bolometric light curve of SN 2003bg compared to the synthetic light curves 
of different models (see text). The inset shows the evolution of the photospheric velocity (as 
shown in Fig. 4) and the prediction of the various light curve models. 
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Fig. 6. — The relation between progenitor mass and synthesised Ni mass in a number of 
Type Ib/c SNe. 
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Fig. 7. — The relation between progenitor mass and explosion kinetic energy in a number 
of Type Ib/c SNe. 



